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Abstract 

The two-dimensional cell-sorting problem is found to be math- 
ematically equivalent to the one-dimensional random walk problem 
with pair creations and annihilations, i.e. the adhesion probabilities 
in the cell-sorting model relate analytically to the expectation values 
in the random walk problem. This is an example demonstrating that 
two completely different biological systems are governed by a common 
mathematical structure. This result is obtained through the equiva- 
lences of these systems with lattice spin models. It is also shown that 
arbitrary generation operators can be written by the spin operators, 
and hence all biological stochastic problems can in principle be ana- 
lyzed utilizing the techniques and knowledge previously obtained in 
the study of lattice spin systems. 
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1 Introduction 



Lattice models are important in biology when one needs to introduce a 
spatial structure. It has been shown that lattice models are often equiv- 
alent to spin models, which are lattice models for magnetic materials. 
Spin models have long been studied as a purely statistical mechanical 
subject, and various models on various lattices with different types of 
interactions have been investigated in detail. Approximation methods 
and techniques for numerical calculations have been extensively devel- 
oped. Some models have been solved and the analytic expressions for 
physical quantities such as energy, specific heat, susceptibility, or cor- 
relation functions have been obtained. Some of the solvable models are 
known to be equivalent to each other, i.e. analytic relations between 
the quantities of the models are found under some specific relations of 
parameters. 

In this study, it is shown that models for biological systems (or 
ecosystems, or organismic systems) are sometimes governed by a com- 
mon mathematical structure, even though the models seem to differ 
completely from other. We show that the two-dimensional cell-sorting 
is equivalent to a generalized one-dimensional random walk in the fol- 
lowing manner. The two-dimensional cell-sorting model on a square 
grid is equivalent to the two-dimensional Ising model which is one of 
the most typical statistical mechanical models for magnetism. Prior 
study have shown the two-dimensional Ising model is equivalent to the 
one-dimensional XY model, the Hamiltonian of which can be inter- 
preted as the generation matrix of the random walk problem with pair 
creations and annihilations. 

In Section 2, the spin models are briefly reviewed. In Section 3, a 
cell-sorting model is reviewed and its equivalence to the Ising model is 
considered. In Section 4, the analytic relations between the adhesion 
probabilities in the cell-sorting problem on a two-dimensional square 
grid and the expectation values in the random walk problem with pair 
creations and pair annihilations are derived. 

In Section 5, relations between the stochastic movements, partic- 
ularly those of the molecular motors, and spin Hamiltonians are con- 
sidered. Molecular motors are protein molecules that are vital to bio- 
logical motion especially for internal material transport such as mus- 
cle contraction, bacterial motion, cell division, intracellular transport 
along the axons of nerve cells, and genomic transcription. Several im- 
portant families of motor proteins such as ribosomes, kinesins, dyneins, 
myosins and RNA polymerase have been identified. They move along 
periodically structured molecular tracks and show stochastic move- 
ments on an one-dimensional template. In this sutudy, I will consider 
the random walk problem with hard-core interaction, the movements 
of ribosomes on a mRNA, and the movements of kinesins. Finally I 
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will show that generation rules for stochastic movements appearing in 
biological problems can always be represented by the Hamiltonians of 
spin models. 



2 Lattice spin models and their equiva- 
lences 

Lattice spin models: Magnetic materials are composed of atoms that 
have their own magnetic moment. The atoms interact with each other 
energetically, and these interactions are completely quantum mechan- 
ical in nature. Spin models are the mathematical models for magnets 
in which spin operators are assigned to every lattice point and interact 
with each other. Quantum spin is a kind of angular momentum and 
defined through the following commutation relations: 

[sf,s«]=isf, [s«,sl}=isf, [sf,sf} = isv, (1) 

where [A, B] = AB — BA, i 2 = —1, and sf , , sf are the x, y, z 
components of the spin operator at site /, respectively. 

Assume the spin at site i and the spin at site j energetically interact 
with each other. Let the interaction energy be — J(sf +s? +As* s*) , 
where J is the coupling constant and A is the anisotropic parameter. 
The Hamiltonian of the system is written as 

# = - J E[ s ^ + s ^ + As ^]> ( 2 ) 

(id) 

which is the sum of the interaction energies over all the spin pairs 
that interact each other. The eigenvalues of H are the energies which 
can be realized in this system. The model ([2} is called the XXZ model, 
its isotropic case A = 1 is called the Heisenberg model. 

Spin operators are classified by the spin magnitude S. Possible 
values of S are 0, 1/2, 1, 3/2 and m/2 where m is a natural number. 
Let us consider the most typical case: S — 1/2. In this case, there 
are two eigenstates of s\ with the eigenvalues 1/2 and —1/2. We call 
them the spin up state and the spin down state \—}i, respectively. 
Generally there exist n = 2S + 1 eigenstates of the spin operator sf 
when its magnitude is equal to 5*. 

Then we find that the energy contribution from the last term in 
([2]) is obtained from the product of the eigenvalues of interacting two 
spins. We call this factor the Ising interaction. Spin model with only 
the Ising interaction is called the Ising model. 

The one-dimensional Ising model 

N 

h = -jY,44+i (3) 
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was first solved by Ising (1925). When we consider the Ising model 
with N sites, there are 2^ possible configurations of spin up and down 
states, because two states are possible for each site i (i = 1, 2, . . . , N). 
The total energy is a function of the configurations. Thus the en- 
ergy can be written as E k (k = 1,2, . . . ,2 N ) where each k denotes 
a configuration. Following the statistical mechanics, the probabil- 
ity to find a configuration with an energy Ek is exp(— Ek/k^T)/Z, 
where T is the temperature, fee is the Boltzmann constant, and Z = 

2 N 

J2k=i ex P(~ Ek/k&T) is the normalization factor called the partition 
function. Then the expectation value of the energy, for example, is 
obtained as 

2 N 

(E) =Y^E k eM-m)/Z, P=t^, 
which can also be written as 



We then introduce the Gibbs free energy F by the relation — f3F = 
log Z. The expectation value of the total energy E, and other quanti- 
ties like magnetization, specific heat and susceptibility, are obtained as 
the derivatives of F. Models are usually said to be solved when their 
Gibbs energies are analytically obtained in a closed form. 

The Gibbs free energy of the one-dimensional Ising model have 
also been obtained in a sophisticated way in introducing the transfer 
matrix method (Kramers and Wannier 1941; Kubo 1943). Let m, be 
the eigenvalue of sf . The partition function of this model is written as 



Z = ^Texp[-/3£ fc ] 
k=i 

= "' exp[/3J(mim 2 + m 2 m 3 + • • • m N mi)] 

= ^ ••• exp(/3JmiTO 2 ) exp(/3Jm 2 m 3 ) • • • exp(/3Jm A rm 1 ), 

where the periodic boundary condition tun+i — mi is assumed. When 
one introduce 2x2 matrix (V) mm i = exp(/3 Jmm'), i.e. (V)i i = 
00-5,-5 = cx P (+/W/4) and (V) i i = (V)_ii = cxp(-/3 J^thc 
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partition function is written as 



Z = E E E (y)n ll m 2 (V) m2 m 3 ---(V) r 

ni— ±-| m 2 — d 

E (Or 



mi- ±i m 2 — ±i m<v-±T 



m=±| 

= TrF Ar = Af + A^ = Af[l + (^) 7V ], 

Ai 

where Ai and A2 are the eigenvalues of the matrix V: \\ = 2 cosh(/3 J/4), 
A2 = 2 sinh(/3 J/4). The free energy per site in the thermodynamic limit 
is 

> = > = ^ log ^ (1 + ( aT )W)] - ^ log Al - 

Effects coming from the boundary vanish when one take the thermody- 
namic limit N — > 00. Hence the free energy is obtained if the maximum 
eigenvalue of the transfer matrix V is obtained. 

When we consider the two-dimensional Ising model on a square 
lattice with the size M x N, we have to introduce a 2 M x 2 M transfer 
matrix. This huge matrix have been diagonalized and the free energy 
have been obtained (Onsager 1944). 

All the spin operators s\ in the Ising Hamiltonian commute with 
each other: [sf,s|] = 0. Then the system is called a classical model 
because there exist no quantum effect coming from non-commutativity 
of operators. 

Next let us introduce sf — sfztisf. From the commutation relation 
(JTJ) , it is derived that sf maps the spin down state to the spin up state, 
s ; ~ maps the spin up state to the spin down state, and otherwise it 
works as the zero operator: 

^ + |->; = l+>;, sT\+)i = \-h =0, *f|-)i = o. 

The sum of the first two terms in ([2]) is equal to (sf sj + s~s^)/2, and 
this term transfers the spin up state from site j to i, or from site i to 

r- 

stsj\-)i\+h = sr s +\+) i \-) j = 1-),!+),. 

Contributions coming from the other states are zero: sfsj\+)i\—)j = 
«^*7l+)i|+)j = ataj\-)i\-) s = and s7s/|-) i |+) j = »r«/l+><l+>j = 
s 7 s j\~)i\~)j = 0- Thus the sum of the first two terms in ((2|) is the 
two-body flip operation, in which + moves from the right to the left, 
or from the left to the right. Therefore the Hamiltonian @ is the sum 
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of the two-body flips and the products of the eigenvalues of interacting 
two spins. 

The spin model interacting via sf sj and sfsj, 

H = - J£[(l + 7)*?«? + (1 - 7)*?^] - ^E s - W 
<*>j> 2 

is called the XY model, where 7 is the anisotropy parameter. The 
factor h in the last term is an external magnetic held applied to the 
z direction, which will become important to consider the equivalence. 
The interaction is a sum of the two-body flips and hence it is not trivial 
to find the eigenstates, in contrast with the fact that the eigenstates of 
the Ising model is obtained immediately In the one-dimensional case, 
the XY model have been solved exactly (Lieb et al. 1961; Katsura 
1962; Niemaijer 1967). 

The Ising model with an external field applied to the x direction, 
— h x J2i s ii i s called the transverse Ising model. The transverse sus- 
ceptibility at h x = for the one-dimensional case with 5 = 1/2 have 
been exactly calculated (Fisher 1960). The transverse susceptibility 
at h x = for the two-dimensional case with 5=1/2 have also been 
obtained (Fisher 1963). The exact free energy for the one-dimensional 
case have been obtained (Katsura 1962; Pfeuty 1970). The transverse 
term is expressed as s x = (sf + s7")/2, and thus it represents inde- 
pendent one-body flip of each spin with the probability proportional 
to h x . The transverse susceptibility at h x = of the one-dimensional 
transverse Ising model with arbitrary spin magnitude 5 have been ex- 
actly obtained (Minami 1996). The susceptibility and the specific heat 
for general Ising type interactions (sf) m (s|)" have also been exactly 
obtained (Minami 1998). 

The Hamiltonians of the XY model and the transverse Ising model 
are the sums of operators which do not always commute with each 
other. Non-commutatibity of operators induce quantum effects, and 
these models are typical examples of the quantum spin system. 

Equivalences: It has been derived that some of these spin models are 
equivalent to each other in the following sense (Suzuki 1971). The two- 
dimensional Ising model is solved using the transfer matrix V, where 
the free energy is obtained from the maximum eigenvalue of V. Let 
|0o) be the eigenstate of V corresponding to the maximum eigenvalue. 
With an appropriate choice of parameters, the Hamiltonian H of the 
one-dimensional XY model commutes with the transfer matrix V of the 
two-dimensional Ising model: [H, V] = HV — VH = 0. Hence these 
models have a common set of eigenstates: H and V can be diagonalized 
simultaneously. In particular, the eigenstate \4>q) corresponding to the 
maximum eigenvalue of V and the eigenstate for the smallest eigenvalue 
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of the one-dimensional XY model are the same. This structure implies 
that the expectation values in the two-dimensional Ising model and the 
expectation values in the lowest energy state (i.e. the ground state) of 
the one-dimensional XY model analytically relate to each other. 

The two-dimensional Ising model and the one-dimensional XY model 
are equivalent to each other, where an anisotropic limit of the XY 
model with an external field is the transverse Ising model. The en- 
ergy of the two-dimensional Ising model is determined from the simple 
product of the eigenvalues of interacting two spins. The interaction is 
simple but the model is not easy to solve because it is defined on the 
two-dimensional square grid. The interaction of the one-dimensional 
XY model is the two-body flip and it is not trivial to find the eigen- 
states, though the lattice is a one-dimensional chain and simpler than 
the square lattice. The interaction of the transverse Ising model is sim- 
ple Ising interaction, however there exists an external magnetic field 
in the x direction, which does not commute with the Ising interaction. 

It is derived that the correlation functions in the two-dimensional 
Ising model and those in the lowest energy state of the one-dimensional 
XY model satisfy the following relation: 

<44c>2D Ising = COsh 2 K^8%)l XY - sinh 2 K*( S V jS V k ) 1D XY . (5) 

The parameters are assumed to satisfy cosh2if* = 1/7, tanh2i4T 2 = 
(1 -7 2 ) 1 / 2 //!, and sinh 2K t smh2K? = 1, Ki = (3J l (i = 1,2), where 
Ji and J2 are the vertical and horizontal coupling constants of the 
square lattice Ising model, respectively (Suzuki 1971, in which the Ising 
model is written in terms of the Pauli operator crj where = crJ/2). 
The expectation value of a quantity Q in the state |</>) is calculated 
using the operator Q which corresponds to Q, and using the expansion 

\4>) = Efe c fel fc ) as 

(Q) = {c/>\Q\<P) = J2 c t'Ck{k'\Q\k), (6) 
fcfc' 

where c* k , is the complex conjugate of Ck> and (k'\ is the dual state of 

W). 

Many other examples of equivalences of lattice spin models have 
been investigated. The six-vertex model, which is a two-dimensional 
lattice model, is equivalent to the one-dimensional XXZ model. The 
eight-vertex model is equivalent to the one-dimensional XYZ model. 
General formula for equivalences between (i-dimensional classical sys- 
tems and (d+l)-dimcnsional quantum systems have also been obtained 
(Suzuki 1976). 
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3 A model for cell-sorting and equivalence 
to the Ising model 

Cell sorting problem: Among the processes involved in the forma- 
tion of an animal, one of the most important phenomena is the self- 
rearrangement of cells leading to the formation of functional structures. 
Starting from a random mixture of cells from different origins, the cells 
reassemble themselves, begin to form clusters of cells of the same type, 
and simulate their normal histological patterns. These movements are 
directed, spontaneous, and proceed in the absence of external forces. 
This pattern formation phenomenon is known as cell-sorting. 

The mechanisms that determine why cells adhere to one another, 
i.e. the forces that drive cell movement during the relevant processes, 
are an important area of research. Steinberg made two assumptions 
to explain certain kinds of cell rearrangements. He assumed that cell- 
sorting required spontaneous progressions of motile and mutually ad- 
hesive cells to configurations that have minimum adhesive free energy 
(Steinberg 1962a; b; c; 1963; 1970) . His assumptions are i) any con- 
tact between cells has an adhesion energy depending on the cell types, 
and ii) cells are mobile and can reach a global energy minimum con- 
figuration independent of their initial condition. These assumptions 
indicate that differential intercellular adhesion and random movement 
of cells are the basic mechanisms of this self-organizing phenomenon. 

This differential adhesion hypothesis has been checked against var- 
ious spatial restrictions, and various additional assumptions on the 
movement of cells have been developed. The cell assumed to move 
on either a two-dimensional square grid or a three-dimensional cube, 
has been investigated analytically, tested through numerical simula- 
tions, and applied to real systems (Goel et al. 1970; Goel and Leith 
1970; Leith and Goel 1971; Goel and Rogers 1978; Rogers and Goel 
1978; Mochizuki et al. 1996; 1998; Mochizuki 2002). Cell movement 
on a hexagonal grid has also been considered (Antonelli et al. 1973; 
1975; Rogers and Sampson 1977). Cells are represented by not only 
hexagons but also general n-gons (Matela and Fletterick 1979; 1980), 
Voronoi polygons (Sulski et al. 1984), or polygonal cells (Graner and 
Sawada 1993). Cells are represented by the large-Q Potts model, in 
which each cell can take Q internal states (Graner and Glazier 1992; 
Glazier and Graner 1993, Nakajima and Ishihara 2011). A viscous liq- 
uid model with interfacial tension was considered (Gordon et al. 1972). 
Steiberg's theory was modified using dynamical equations of a molec- 
ular nature, and cell-sorting was found to occur in a near-liquid state 
(Greenspan 1981). A continuous mathematical model was proposed to 
analyze cell-sorting in Dictyostelium discoideum (Umeda 1989; Umeda 
and Inouye 1999; 2004). 
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In the present study, I concentrate on the model introduced by 
Mochizuki ct al. (1996) on the two-dimensional square lattice. This 
model is directly equivalent to the Ising model, which has been exactly 
analyzed as a model for magnetic materials. 

Model and the equivalence: Let us consider two kinds of cells dis- 
tinguished by color: black and white. The cells are assumed to form a 
regular square lattice. Let Abb be the strength of adhesion per cell con- 
tact between black and black cells, and Aww an d Abw be the strength 
between white and white, and black and white cells, respectively. One 
can estimate the total adhesion A& when the configuration k of the 
black and white cells is known: A& is the sum of all the strength of 
adhesion between cells. Cells exchange their positions between nearest 
neighbors. Let m be the rate of exchange of the location. It is assumed 
the probability of finding a configuration having total strength equal 
to Afe is proportional to exp(Afc/m). 

Let us introduce the differential adhesion A = Abb + Aww — 2 Abw- 
The total adhesion A& changes its value by A as an unit when two 
cells change their locations, and hence the system is controlled not 
through the energies Abb , Aww an d Abw independently, but through 
the differential adhesion A. Typical equilibrium configurations are 
calculated by Mochizuki et al., where one can find a pattern in which 
black and white cells appear situated like a checker-board for A/m = 
—2, the pattern seems almost random when A/m = 0, and it appears 
segregated when A/m is positive and large. 

This model of cell-sorting is equivalent to the Ising model when 
Abb = Aww- One can introduce a direct correspondence in which a 
black cell at site i corresponds to the spin up state \+}i, and a white 
cell at site i corresponds to the spin down state |— )j. The interaction 
energy between two up states, or two down states, is equal to — J/ 4, and 
that between up and down states is equal to +J/4. Let us introduce 
a constant Jo and set 

Abb = Aww = — (— J/4+J ), 
Abw = -(+J/4 + J ). 

The overall minus signs are introduced because the probability to find 
a configuration with total adhesion strength A& is proportional to 
exp(A/j/m) in the case of the cell-sorting model, though the proba- 
bility to find a configuration with total energy is proportional to 
cxp(— Ek/k^T) in the case of the Ising model. It is easy to check that 
the contributions from J cancel in the calculations of expectation val- 
ues. The Hamiltonian is H = — J^(ij) s i s j + \zNJq, where TV is the 
total number of sites and z is the number of interacting pairs (ij) from 
each site i, e.g. z = 2 in the one-dimensional linear chain and z = 4 in 
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the two-dimensional square lattice. The parameters in the cell-sorting 
model and those in the Ising model arc related by A — J and m — k-gT, 
and hence A/m = J/k^T. The probability of re-arrangement increase 
as the temperature increase. The correspondence is consistent with 
the facts that J is the coupling constant which determines the unit 
of energy scale in the Ising model, and A is the the difference adhe- 
sion which determines the unit of the total adhesion in the cell-sorting 
problem. 

As a model of magnetism, A/m = —2 indicates that the Ising model 
is an antiferromagnet. The checkerboard- like pattern is that called the 
Neel order. The case with A/m = corresponds to a paramagnet, 
where there is no interaction between spins, and each spin flips inde- 
pendently. The cases with positive A/m correspond to ferromagnets, 
where the spin of the nearest-neighbors tend to become parallel and 
tend to form clusters. In the case studied by Mochizuki et al., the 
number of the black cells and the number of the white cells are fixed 
to be the same. 

The two-dimensional Ising model has a critical point T c at a non- 
zero and finite value of temperature. Below T c , clusters of spin up 
states appear, which corresponds to the existence of cell-sorting phe- 
nomena, whereas above T c there is no cell-sorting. In the case of the 
one-dimensional Ising model, the critical temperature is equal to zero. 

The quantities studied in the cell-sorting problem correspond to 
quantities in the Ising model. Let us assume the periodic boundary 
condition. Following the notations in Mochizuki ct al., let pe be the 
fraction of the black cells and pw be the fraction of the white cells, then 
Pb + Pw = 1 • Let (jbb be the fraction of black cells in the neighborhood 
of a randomly chosen black cell, <7ww the fraction of white cells in the 
neighborhood of a randomly chosen white cell, qwB the fraction of 
white cells in the neighborhood of a randomly chosen black cell, and 
<Zbw the fraction of black cells in the neighborhood of a randomly 
chosen white cell. Then gsB + 9wb = 1 and gww + <7bw = 1 • The two- 
body correlation function in the Ising model is the expectation value 
of the product of eigenvalues ±1/2 of corresponding two sites. Thus, 
by definition, the nearest-neighbor correlation function is written using 
the probabilities as 

11 1-1 -11 -1-1 

\ s i s j/Ising — ^PB QBB + -^-^-PB QWB + —-^PW 1BW + ——PW qwW 
1 , N 1 

= ^Pb (2<?bb - 1) + -^Pw (2g W w - 1) 

= ^( 2 Pb <?bb + 2p w 9ww - !)• (7) 

The magnetization (sf ) ising is the expectation value of the eigenvalue 
associated to each site, and obviously is equal to \pb + {—^)pw- How- 
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ever, it can also be written in terms of the probabilities pb <Zbb and 
Pw <7ww- The probability that the two cells in a randomly chosen 
nearest- neighbor pair are both black is pb ?bb > the probability that one 
of the two cells is black and the other cell is white is pb qwB + Pw Qbw > 
and the probability that the two cells are both white is pw <Zww • Con- 
sidering the average magnetization per pair, we obtain 

(sf)lsing = ^[(^ + ^PB 9BB + (^ ~ 7j)(PB <?WB +PW 9Bw) 

+( _ 2 ~~ 2^ W qww ^ 

= ^(PB<?BB-PW9WW) (8) 

From ([7]) and (|8]l. pb<7bb and pw9ww are expressed by the correlation 
functions in the Ising model as 

PB <?BB = (sf S|)lsing + J + («f )lsing, 
PW ?WW = (sf s |)lsing + J" (sf)lsing! (9) 

This argument can be generalized to obtain probabilities for other 
types of spin pairs. When the site i and j locate in the next-nearest- 
neighbors (or in other specific locations) each other, ^bb etc. should 
be redefined as the fraction to find cells of corresponding colors in the 
next-nearest-neighbors (or in the other specific locations). Then (J9j) is 
still valid. With the open boundary, © is also valid except boundary 
terms, which vanish when one assume the periodic boundary condition. 

In the one-dimensional Ising model, there is no spontaneous sym- 
metry breaking, and hence (sf)i s in g = 0, thus pb = 1/2. The nearest- 
neighbor correlation function of this case is immediately calculated by 
the transfer matrix method and one obtains qBB = 2(tanh(J/4A;T)/4 + 
1/4 + 0) = (tanh(A/4m) + l)/2, which is consistent with the result ob- 
tained for the cell-sorting model. 



4 Equivalences between biological systems 

The cell-sorting model with adhesion strength is equivalent to the Ising 
model in each dimension. The two-dimensional Ising model on the 
square lattice is equivalent to the one-dimensional XY model with an 
external magnetic field, and it is also equivalent to the one-dimensional 
transverse Ising model. The XY model and the transverse Ising model 
are equivalent to systems with two-body flip, or independent one-body 
flip, with the energy estimated as the product of eigenvalues of two 
states. Next I consider a stochastic process that corresponds to the 
XY model. 
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Let us consider identical particles on the one-dimensional lattice. 
Each particle locate on a lattice point x, where x is integer. The 
particles move right or left in each step with the rate pr and pl, re- 
spectively. The particles interact with the hard-core interactions: they 
cannot move into the space already occupied by other particles. Let 
us introduce the following notation 



/j+i 



etc. 



Operators sj and s^ +1 work nontrivially only on the j-th and j + 1-th 
sites as 



1 1 1 1 



(10) 



When we regard + as a particle and — as an empty site, operations 
(|10[) are the mappings in which a particle moves to one of its nearest 
neighbors. The matrix representations of the operators are 
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(11) 



where the basis set used to represent the matrix is 



{|+>il+>*+i. \-h\-)j+i, \+)j\-h+u l->il+)i+i,}- 

The matrix [ operates on the site j and j + 1, and operates as an 
identity on the other sites. Then the contributions from the random 
hoppings arc written as 

N 
N 

- E 

Let {x\, X2, ■ ■ ■ , x mp } be a configuration of the particles, in which the 
sites xi, X2, ■ ■ ■ , x m are occupied and the other sites are empty. Let 
{xi]k (k — 1, . • • , 2 N ) be the possible 2 N configurations, where each 
k denotes one configuration. Let P({xi}k;n) be the probability that 
the configuration of the particles is {xi}k after the n-th step. Let us 
introduce the probability vector P n , in which the /c-th element of P n 







pl 

PR 



(13) 



ii'4-1 
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is (P„)/t = P({xi}k]n). Then the random hoppings are generated by 
P„+i = HP n . This process conserves the number of particles. In the 
case of pr, = pl = 1/2, the operator (fT2"|) is reduced to the Hamiltonian 
of the XY-model Q with 7 = and h = 0. 

When one introduce nearest-neighbor pair creation and annihila- 
tion processes with the probability pu and pd, respectively, and also 
multiply the weights —h/2 and h/2 for the nearest-neighbor pair par- 
ticles and pair of the empty sites, respectively, H becomes 



H 



N 

E 

3=1 



-h/2 

Pb 





h/2 









PR 






PL 





This operator (|14p is written by the spin operators as 



JV 



H 



3 



PLSj S j+1 



-Pus, s j+1 



PDSj S j+1 



3 = 1 



(14) 



+ *?+!)]■ 



When pr = pl = 1/2 and pu = Pd = 7/2, (|T5|) is reduced to the 
Hamiltonian of the anisotropic XY-model (j4j with an external field h. 

The first 2x2 block element in (|14|) represents the following pro- 
cesses. Two particles meet and annihilate with the rate pd, a nearest- 
neighbor pair of particles are created with the rate pu- The weight 
of each state varies spontaneously: multiplied by the factor which is 
the sum of —h/2 and h/2 associated to each nearest-neighbor pair of 
particles and each nearest-neighbor pair of empty sites, respectively. 

Let pu — pd =7/2, and assume 7 > 0. All the matrix elements in 
the first block in (|14p can become non-negative by a rotation as 



R- 



-h 
7 



R = X 



33 + 1 



33 + 1 



where A = ^7 2 + h 2 /2 and 



R 



cost 
sin f 



— sm( 
cos 9 



33+1 



tan 9 



/1-7 + ^7 + h 2 n 
h + 7 + V7 2 + ft, 2 4 



Thus the pair creation and annihilation processes with the sponta- 
neous change —h/2 and h/2 can also be regarded as a simple reflection 
multiplied by the factor A, in the above two rotated basis. 

Therefore the Hamiltonian of the XY-model can be regarded as the 
generator of the following two processes, 1) random walk process with 
hard-core interactions, and 2) pair creations and pair annihilations 
with spontaneous change of weights which are determined from the 
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number of pair particles and "pair vacuums". The second process 2) 
is also regarded as the rotated reflections with the rate y/-f 2 + h 2 /2. 

In the relation ([5]), the expectation values ( )idxy are estimated 
in the ground state \<j>o) of the XY model with the coupling constant 
J > 0, and thus — J < 0. The state \4>o) is the eigenstate corresponding 
to the maximum eigenvalue of f|15[) because the sign of the coupling 
is now — J = 1 > 0. When we consider the matrix H + cl where 
I is the unit matrix and c > is sufficiently large, all the matrix 
elements of H + cl are non-negative, and the eigenstates of H + cl 
are simultaneously the eigenstates of H. From the Perron-Frobenius 
theorem, it is derived that all the coefficients in \4> ) are non-negative. 

Let Ao be the maximum eigenvalue of H. The maximum eigenvalue 
of the matrix H ~ H/Xo is equal to 1. Let A; be an eigenvalue of H and 
| (pi) be an eigenstate corresponding to A^: H\<j)i) = Aj|0j). Then all the 
eigenvalues satisfy |A;| < |A | = 1. Therefore the state |0o) survives in 
the limit where the operator H is iteratively applied: /f"|0o) = |</>o)- 

In the one-dimensional XY model, there exist a critical field h c , 
which corresponds to the critical temperature T c of the two-dimensional 
Ising model. The ground state is unique for h > h c and two-fold de- 
generate for h < h c in the thermodynamic limit N — > oo. In the case of 
the two-dimensional Ising model, the eigenstate of the transfer matrix 
corresponding to the maximum eigenvalue is unique for T > T c and 
two- fold degenerate for T <T C , and there exists a spontaneous symme- 
try breaking and the magnetization becomes non-zero below T c . This 
implies that the global minimum of the free energy is realized with 
Pb = 1/2 for T > T c , and with p B ^ 1/2 for T < T c . The correspond- 
ing cell-sorting model below T c is, therefore, a generalized model which 
includes the effect of gathering of one kind of cells from outside. 

If the process represented by H is fully stochastic, the matrix H 
should satisfy the conservation of probability: (H)ij > and J2i(H)ij = 
1. The matrix H does not satisfy this condition. However, the state 
|^o) is a steady state and the conservation of probability is satisfied in 
\cj) ). Let 

|<fo ) = ^2 Pq{xi,x 2 , ■ ■ . ,x mp ;n)\x 1 ,x 2 , . . -,x mp ), 

where \x\,X2, ■ ■ ■ ,^m p ) is the direct product of the states {j — 

x 1 ,X2,---,x mp )and \-)j (J j£ x 1 ,x 2 ,...,x mp ), i.e. the state where 
the sites x±, . . . , x mp are occupied and the other sites are empty. The 
coefficient Pq(xi, x 2 , ■ ■ ■ , x mp ; n) is the probability that the system is in 
the state \x%, x 2 , ■ ■ ■ , x„ lp ). The probability distribution {Pq(xi, x 2 , . . . , x mp ; n)} 
is invariant under the operations of H and satisfies the conservation 
of probability. This means that H is stochastic in the subspace corre- 
sponding to the maximum eigenvalue. 
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In the case of the cell-sorting problem on the square grid, from 
(©, ©, and from the relations (sf)2Dismg = (1/2)pb + (-l/2)pw> 
Pb + Pw = 1, and cosh 2^ = 1/7, one obtains 

PB<7BB = (sfjSf fc )2DIsing + T + (sf )2DIsing 

= cosh 2 if* (sp%) id xy - sinh a ifr<s^)iDxy + ~ + (pb - ~) 

= PB + ^[(4 S fe + S J S k )lDXY + ^(4 S fc + S J S fe)lDXY - 1] 

(16) 

The probability <7bb m f| is that for the two cells locating 
and (i,k). The parameter 7 is specified from the restriction K\ — 
K 2 = f3J = A/m. The expectation values are calculated in principle 
from ( • )idxy = (0o I • \4>o), the expansion ©, and the normalization 
conditions (+|+) = 1, (-|-) = 1, (+|-) = and (-|+) = for each 
site i. The magnetization (sf)2Disin g below T c have been calculated 
by Yang (1952). The correlation functions (s^s^) 1DX y {I = x,y) have 
been analytically calculated in the context of the one-dimensional XY 
model (McCoy et al., 1971; Tonegawa, 1981), and (sQ s^ fe )2D ising have 
also been calculated in the context of the equivalent two-dimensional 
Ising model (Wu et al. 1976; McCoy et al. 1977), both for arbitraly j 
and k. In conclusion, the adhesion probability in the two-dimensional 
cell-sorting problem is analytically expressed by the expectation values 
in the steady state of the one-dimensional generalized random walk 
problem. 



5 Correspondences of stochastic processes 
with lattice spin models 

The matrix H (and H) corresponding to the XY model is not itself 
fully stochastic for all the possible states of the particles. In this sec- 
tion, first I formulate the most general form of the stochastic matrix 
for the one-dimensional random walk process with creations and anni- 
hilations, and next consider three examples to show the close relation 
between biological systems and lattice spin systems, and finally prove 
that the stochastic processes in biological systems can always be writ- 
ten in terms of the spin Hamiltonians. 
Let us consider the operator 

JV 

H = J2 



w 3 p v 

Pb W4 

Wl PL 

pr w 2 



(17) 



15 



The element j»d in the matrix [ is the rate of the pair annihilation 
^ I - The element W3 is the rate of the pair to be 
invariant H> an d the other elements are defined 

as the rates of the following transitions: 

Pu ■ ^ w i ■ ^ |->i|-)j+i 

Pn ■ w x : |+)i|->j+i h> 

First let us assume the periodic boundary condition, and consider the 
case Wi — (i — 1, 2, 3, 4). Starting from the state |+)i|+)2 • ■ ■ \+)n — 
! + +■■•+), N kind of pair annihilations which result in the states 

with successive two minus, | hH h), | H 1 h), etc., 

are possible. Thus the conservation of probability is satisfied with 

the condition Npn = 1. Starting from other states such as | H h 

+ •••+), I H (■■■•+), H 1 — ■••+}) etc., it is straightforward 

to convince that two pd's (or two pjj's) are always replaced by a pair 
formed by j?r and pl in the calculations of the total probability, and 
thus the conservation of probability is satisfied with 2po = PR + Pl, 
2fu = Pn + Pl- Similarly in the case of pd, the condition for pu is 
Npu = 1. Thus we obtain p D = pu — 1/N, and p R + Pl = 2/JV. Let 
us introduce the rates u>.; (i = 1, 2, 3, 4), then p-Q, pu, Pr, and Pl should 
be replaced by pr> + W3, Pu + Wi, PR + wi, and pl + W2, respectively, 
in the calculations of the total probability. Therefore the conservation 
of probability is satisfied iff 

1 , , 

Pn+w 3 =pu + Wi = (18) 

{pR + w 1 ) + { PL +w 2 ) = A (19) 

When we assume the open boundary condition, s z N+1 ^ sf, then the 
number of spin pairs (i,i + l) is iV, and the condition (|18p remains true 
though the number of sites now being N+l. The condition (|19|) should 
be satisfied for 2 < j < JV — 1. At (j, j + 1) = (1, 2) and (N, N + 1), 
the probability is conserved iff 

Pr + W! =pl + w 2 = —. (20) 

The condition (fl9|) is satisfied when we assume (|20|). Thus fp~8|) and 
((2U)) are the conditions in the case of the lattice with open boundary. 
These are the general form of the one-dimensional stochastic matrix 
for random walk processes, with pair creations and pair annihilations. 
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The matrix (|17|) is written by the spin operators as 

N 

H = ^2[prsJs+ +1 +PLsjsj +1 +pvsfsj +1 +PdsJsJ +v 

3 = 1 

+As* jS * j+1 - ^ + s" j+l ) + col + C -{s] - 8' j+1 % (21) 
where / is the identity operator and 

A = 1U3 + W4 — wi — W2, h — W4 — W3 

c = i(w 3 + w A + w 1 + w 2 ), c = w 1 -w 2 - (22) 

When one assume the periodic boundary condition s^ +1 — sf, 

the last term in (f2Tj) vanishes: X)jLi( s ) — s i+i) = 0, and thus H is 
independent of c. It is easy to check that 103, and w\ + W2 do not 
depend on c, and hence the process remains stochastic, i.e. satisfy (ITS)) 
and (1191) . with c being a free parameter. This is a result which comes 
from the translational invariance of the system. 

Next let us consider one-body creations and annihilations of the 
particles. Let us introduce an operator 

nj = ~+s*. (23) 

In our notation, the eigenstate of s| with the eigenvalue +1/2 (the 
eigenvalue —1/2) corresponds to the state with a particle (without a 
particle) at site j. The operator rij takes the value 1 or when the site 
j is occupied by a particle or the site j is empty, respectively. Thus 
rij is called the number operator. Let us consider the case where a 
particle is added at site j. The new state with the added particle is 
created by s^~, and the old state without the particle is removed by 
-(1-nj): 



" 


1 " 




' 


" 








j 





1 



(Note that these operators work non-trivially only in the subspace 
where rij =0, states in which the site j is empty, and otherwise they 
work as the zero operator.) Thus the term AHj' = aj[s^ — (1 — rij)] 
should be added to the Hamiltonian H in (|2"Tj) . where aj is the rate of 
creation at site j. Similarly, if a particle at site j is annihilated with 
the rate f3j, the term AiJ~ = /3j[sJ — nj] should be added to (|2"Tj) . 
(These operators work only in the subspace where rij = 1, states with 
a particle at site j.) Because = sj ± is|, these terms are the trans- 
verse and parallel magnetic fields applied to the site j. Thus we find 
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that external fields induce the spontaneous creations and annihilations 
(attachments and displacements) of the particles. 

Random walk with hard-core interaction: From (flU)) and (fT5|) 
- ([2"2")) . the fully stochastic random walk process with hard-core interac- 
tions is written as 



where rfa/N = pr, rj^/N = j»l and 7/r + t/l = 2. This operator is 
written by the spin operators as 



When rjR = rjh = 1, the operator (f25|) is the Hamiltonian of the XXZ 
model ([2|) with A = 1, i.e. the Hamiltonian of the Heisenberg model. 
This case is a kind of the hard-core boson system, which is originally 
proposed as a model for helium superfluidity (Matsubara and Matsuda, 



The anisotropy rj^ ^ rj^ can be removed by an unitary transforma- 
tion (Henkel and Schutz 1994, in which the eigenstate with the eigen- 
value — 1/2 is regarded as the state with a prticle) provided jtri/l 7^ 0. 
Let 




10 

10 

r] h 

m 



(24) 



N 



1956). 



N 




(26) 



where rij is given by (|23|) . The commutation relation [rij, s^] = 
and the expansion 




yield Vsfy- 1 



q ±rj 'Sj and VsjV 1 = s|. Let 5 = y/rjij rjR and one 



obtains 



-1 



D^fo" + ?4*7 + t) + + \i\v- 1 
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Hence it is derived that the anisotropy of the hopping rates in the ran- 
dom walk process can be handled as the anisotropy A = 1/^/?7r?7l of 
the quantum coupling in an uniform spin chain. 

Ribosome moving on mRNA: A ribosome is a large, complex 
molecule that synthesizes a protein molecule using the genetic mes- 
sage coded on mRNA as the template. RNA comprises four kinds 
of nucleotides, and a triplet of nucleotides constitutes a codon. Each 
possible type of codon corresponds to one species of amino acid: 61 
kinds of codon lead to 20 species of amino acid, whereas three special 
codons indicate termination of translation. The information enclosed 
in the codon sequence is translated by the ribosome into the amino 
acid sequence of the encoded proteins. 

A ribosome binds to an mRNA and begins to synthesize the pro- 
tein by adding an amino acid (this is referred to as initiation). After 
a biochemical reaction for the elongation of the protein, the ribosomes 
moves forward on the track by one codon (i.e. one elongation has oc- 
cured) . Finally, the ribosome reaches the termination codon and leaves 
the mRNA, releasing the protein (this is referred to as termination). 

A number of ribosomes can be simultaneously attached to one 
mRNA template. A ribosome can move forward on the track pro- 
vided that the next codon is not captured by another ribosome, i.e. 
ribosomes are moving on the mRNA template interacting via hard- 
core interactions. A ribosome is often treated as a molecular motor, 
and this collective movement process along the mRNA chain shows 
a correspondence with a one-dimensional driven lattice gas, or with 
vehicular traffic on a road. 

MacDonald et al. introduced a stochastic process now known as the 
asymmetric simple exclusion process (ASEP), as a model for the move- 
ment of ribosomes on a mRNA (MacDonald et al. 1968; MacDonald 
and Gibbs 1969). The model was first introduced in the biophysical 
literature, and later studied from a purely theoretical viewpoint (see 
for example Derrida 1998). 

The actual movement of a ribosome is closely coupled to its internal 
mechanochemical processes to synthesize a protein. Accounting for 
these processes, the ASEP has been generalized to have seven (Basu 
and Chowdhury 2007), five (Garai et al. 2009), or two (Ciandrini et 
al.) distinct biochemical states in each cycle. The ribosome movement 
is also characterized by a pause and translocation, which defines the 
time of its residing at a corresponding codon. The ASEP has been 
generalized to have one or more slow codon bottlenecks (Kolomeisky 
1998; Chou and Lakatos 2004; Dong et al. 2007a; Dong et al. 2007b). 
The ASEP with different hopping rates associated with each site has 
been considered (Shaw et al. 2003; Shaw et al. 2004; Romano et al. 
2009). A ribosome recycling mechanism has been introduced (Chou 
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2003), in which a part of the ribosome detaches at the termination 
site and part of them diffuses back to the initiation site. The ASEP 
comprising open boundaries with random particle attachments and 
detachments was also introduced (Parmeggiani et al. 2003, Pierobon 
et al. 2006). A stochastic model with a secondary structure of mRNA 
was introduced by von Heijine et al. (1977). 

Now, let us introduce the ASEP, which is a random walk process 
with specified moving rates and with a continuous time variable. Let 
us consider an one-dimensional lattice. Each site j (J £ Z) is empty or 
occupied by a particle which corresponds to a ribosome. Each particle 
stochastically move forward or backward interacting via the hard-core 
interactions i.e. a particle at site j moves to the site j + 1 with the rate 

provided that the site j + 1 is empty, and a particle at site j + 1 
moves to the site j with the rate p\, provided that the site j is empty. 

Let us first assume the periodic boundary condition and consider 
the stochastic matrix (TT7|) with pu = pd = 0, and thus from ([15]) 
W3 = W4 = l/N. Because H is independent of c, the rates w\ and W2 
can be taken as 



PL, 



(c = p h -p K - —). 



These rates satisfy (fT8|) and (|19[) . The corresponding stochastic matrix 
is 



Has 



EP — 



1 N 

-Y 



10 

10 

(l-5)- m , VL 

7 7R (1 + S)- Vl 

1 + AHasep, 

iv 



j'i+i 



AH 



ASEP 



2 

N 



1, _ 1 



3=1 



i)P7) 



This is the ASEP with discretized time step. Because of the peri- 
odic boundary condition, the last term vanishes and the operator is 
independent of 5. The rates 6 and —6 always appear pairwise in the 
summation X)j=ij an d finally give no contribution to the probability. 
When one set ijr = 1 — S and t]l = 1 + 5, the system is reduced to 
and (1231) . In this case, the ASEP is a simple random walk process 
with hard-core interactions. 

As a model of the movement of ribosomes, the rates are asymmetric 
and the open boundary condition should be introduced. With these 
restrictions, the expression of 7?asep by spin operators have already 
been written by Sandow (1994, see also Alcaraz 1994). The condition 
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(|20| is satisfied provided S = 0. The attachment at the site j = 1 and 
the displacement at the site j = N + 1 are introduced by the terms 
Ofi[si — (1 — rii)] and /3jv+i[s iv+1 — n/v+i]- Then with the use of the 
transformation V = exp[(log q) ^2^=1 3 n i\i which is slightly modefied 
from (f26|) . one obtains 

2 N 1 1 

i=i 

+ai[gsf - (1 - m)] + /3Ar+i[g _(Ar+1) s w+1 - njv+i]. 

Therefore the ASEP with asymmetric rates and with the open bound- 
ary condition is nothing but the XXZ spin chain with A = (q + q~ 1 )/2 
with boundary magnetic fields. Here we find that q is the parameter 
in the quantum group symmetry of the XXZ spin chain (Pasquier and 
Saleur 1990; Jimbo and Miwa 1993). 

Other generalizations of the ASEP reviewed above can be handled 
as generalizations of spin Hamiltonians. The biochemical states can be 
introduced using spin operators with the spin magnitude S, in which 
the number of the eigenstates is n = 2S + 1, and the bottleneck or 
the non-uniform hopping rates correspond to non-uniform coupling 
constants of spin chains. The ASEP with continuous time variation is 
also governed by the same operator (|27l) : its time dependence is 

jV{t) = (AF ASEP )P(t), 

where each element of P(t) is equal to P(x\ ,x%,..., x mp ; t) , the proba- 
bility that the system is in the configuration {x\, X2, ■ ■ ■ , x mp } at time t. 

Kinesins moving on a microtube: Kinesin is a single molecular 
motor observed in vitro to move along a linear microtubule template. 
It moves stochastically and stepwise along a one-dimensional track (for 
example, a review by Yildiz and Selvin 2005). 

Experimental methods to measure the biochemical and biomechan- 
ical properties of a single-molecule enable us to observe the movement 
of a single kinesin. It was shown that kinesin moves stepwise along mi- 
crotubules occasionally moving both forwards and backwards (Kojima 
et al. 1997). It was observed that kinesin is released spontaneously 
from the microtubule (Block et al. 1990), and that the increase of 
the load results in an increasing rate of dissociation (Coppin et al. 
1997). It was suggested that two or more sequential processes dom- 
inate the biochemical cycle (Svoboda et al. 1994), and the time to 
force generation after release of ATP was measured (Higuchi et al. 
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1997). Force- velocity curves was obtained for single kinesin molecules 
(Svoboda and Block 1994). 

Pioneering theoretical models related to the dynamics of kinesin 
are already known. An elementary "barometric" relation was intro- 
duced for the driving force (Fisher and Kolomcisky 1999a; b). Nearest- 
neighbor kinetic hopping models with arbitrary forward and back- 
ward periodic rate constants with three generalizations was introduced 
(Kolomeisky and Fisher 2000a). The one-dimensional random walk 
process with general waiting-time distributions, finite side branches, 
and annihilation was also considered (Kolomeisky and Fisher 2000b). 
The observed movement of kinesin can be described adequately by 
simple discrete-state stochastic models (Fisher and Kolomeisky 2001). 
The model has also been applied to analyze the dynamics of myosin- 
V (Kolomeisky and Fisher 2003). Brownian particles moving in two 
or more periodic but spatially asymmetric and stochastically switched 
potentials have also been considered (Jiilicher et al. 1997) The models 
have been summarized in review articles (Reimann 2002; Kolomeisky 
and Fisher 2007). 

Let us consider a model for kinesin on a microtubel (Kolomeisky 
and Fisher 2000a). A kinesin is assumed to locate on a site j of a 
one-dimensional lattice, where j = 1, 2, . . . , N + 1. Let us introduce an 
index k (k = 0, 1, 2, . . . , K — 1) which distinguishes the internal states 
of kinesin. A kinesin changes its internal state from k to k + 1 with 
the rate Uk, and from k to k — 1 with the rate Wk- A kinesin at site 
j in the maximum internal state k = K — 1 can move to the next site 
j + 1 and initialize its state as k — with the rate uk-i- A kinesin at 
site j + 1 in the minimum state k = can move to the previous site j 
and have the maximum state k = K — 1 with the rate wq. 

Let us consider the case with two internal states: the case K = 2. 
Each site takes one of the three possible states: k = 0,1, and "empty" . 
Thus let us introduce the spin operator with the spin magnitude 5=1, 
in which there are 25+1 = 3 possible eigenstates with the eigenvalues 
S z = —1, +1, and 0. Let us assume that the states with S z = —1, +1, 
and correspond to the states k = 0, 1, and "empty", respectively. 
Let the generating matrix be Hkf = I + AHkf, and consider AHkf- 
The hopping from the state k = to k = 1 at site j is expressed as 
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The backward hopping from the state k = 1 to is, in the same way, 
expressed as HJ — wi(sjsj — Ij). The hopping from the state k = 1 
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at site j to the state k = at site j + 1 is expressed as 



Ui 



= Ui 
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L ( 7f. s D( s I+i^=^+i) - (4 - ^ s 3 7f4)(^+i - (4+i) 2 ) 



i 



= -^i [(^ s ?)( s I+i s 7+i) + « S D 2 + (4+i) 2 )] 



where we made use of the relations (sj) 2 + (s|) 2 + (s|) 2 = 5(5 + 1) = 2 
and 



V2 3 V2 



The backward hopping from the state k — at site j + 1 to the state 
= 1 at site j is, in the same way, expressed as 

= -^o + & ~ (-5) a )((«5+i) a - • 

These operators constitute AHkf ■ The operator Hkf is then obtained 
as 

AT+l JV 

i/KF = / + E + + £(^+i + H ij + i)- 

3=1 J=l 



The operator should be normalized by the maximum eigenvalue Aq. 



General correspondence: In these three examples the generation 
of the biological stochastic systems are written in terms of the spin 
Hamiltonians. The correspondence between stochastic systems and 
spin models are more general. Let us consider a system with n discrete 
internal states. The law governing changes from one state to another 
is written as a matrix of order n. Let us derive the fact that this 
matrix can always be expressed by the spin operators with spin 5, 
where n — 25 + 1. 

Let us consider the unit matrix I and the matrices (s z ) k (k € N). 
They arc diagonal and the diagonal elements of (s z ) k are 5 , (5 — 1) , 
. . . , (-(5-l)) fc and (-S) k . The set of the matrices I, s z , (s z ) 2 , (s z ) 2S , (s z ) 2S+1 
are not independent because they satisfy the eigenequation of s z , which 
is a polynomial of order 25+1. For the purpose to show the indepen- 
dence of /, s z , (s z ) 2 , . . . , (s z ) 2S , let us consider the following determi- 
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nant 



(IS S 2 
1 S-l (S-l) 2 



S 2S 
(S-l) 



det 



1 



1 -(S-l) (-(s-i)) 2 

\ 1 -S (-S) 2 



(-(s-i)) 



J 



l<i<j<n 



where n = 25 + 1 and 



Z\ = S, Z2 = S — 1, . 



-S 1 . 



This is the Vandermonde's determinant and, in our case, clearly non- 
zero. Hence all the diagonal matrices can be expressed as a linear 
combination of I and (s z ) k (k = 1, 2, . . . , 2S + 1). 

In the case of spin operators with the spin magnitude 5", it is derived 
from the commutation relations that the eigenstate of s z corresponding 
to the eigenvalue M satisfies s ± \M) = [S(S+1)-M(M±1)] 1 / 2 \M±1). 
Thus the matrix elements of satisfy (s ± )ij 7^ if and only if j = i±l, 
and ((s^) l )ij 7^ if and only if j = i ± I. Let us introduce a diagonal 
matrix P k by (Pk)ij = 5 ki S kj , i.e. (Pk)kk = I and (Pk)ij = when 
i 7^ k or j 7^ k. It is easy to convince that 



and all the other matrix elements are equal to 0. Therefore arbitrary 
n x n matrices can be expressed in terms of the spin operators with 
spin S, where n — 2S + 1. 

The Hamiltonian generally has the form 



where a and a' denotes z, +, or — . An external field can be introduced 
through the terms with the powers (m, 11) = (0, 1) or (1, 0). Therefore 
all the biological systems with finite number of states and discrete law 
of change can find its equivalent spin model. 

6 Summary and discussions 

In this study, it was derived that the adhesion probabilities in a two- 
dimensional cell-sorting model are analytically expressed by the ex- 
pectation values in the one-dimensional random walk model with pair 



(Pk(s+)% ^ i = k, j = k + l, 
(Pk(s-) l h ^ i = k, j = k-l, 




(28) 
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creations and annihilations. I believe that this is the first example 
in which completely different biological systems show mathematical 
equivalence to each other. It was also derived that the equivalencies 
with spin systems are general, i.e. the generation rules of stochastic 
movements in biological systems can always be written in terms of the 
lattice spin Hamiltonians. This provides a path to study biological 
systems using the techniques and results already obtained in the area 
of spin systems. 

It should be noted that general correspondences between d-dimensional 
quantum spin systems and (d + l)-dimensional Ising-type spin systems 
are known to exist (Suzuki 1976). This implies that for each stochas- 
tic system in cZ-dimension with hopping, creation and annihilation of 
particles, there exist equivalent (d + l)-dimensional model with an 
adhesion- type structure. 

Although the stochastic models can be expressed in terms of spin 
operators, the corresponding Hamiltonians are usually not simple. In 
such situations, numerical calculation techniques might be powerful, 
and many sophisticated techniques have been developed for the lattice 
spin models (see for example the series by Domb and Green 1972) . 
Exact results also exist, for example, in the case of a = a' = z in the 
one-dimensional case of (|2"8|) . that are n-state Ising-type systems with 
independent one-particle creations and annihilations (Minami 1998). 

The dynamics of the Ising model have been studied through the 
system proposed by Glauber (Glauber 1963) or Kawasaki (Kawasaki 
1966). The cell-sorting model considered in this study with pair ex- 
change probability corresponds to the Kawasaki dynamics. The dy- 
namics of n-spin flip models were considered, in which the two-spin 
flip-model is equivalent to the XY chain apart from the boundary 
terms (Felderhof and Suzuki 1971). It is interesting to consider the 
equivalencies of dynamical properties. 

The XXZ model is related to the random walk model, and there- 
fore related to the Brownian motion. The one-dimensional XXZ model 
is equivalent to the six-vertex model, and it is known to demonstrate 
fractal structure in its configuration space (Minami 2010). The XY 
model with periodically varying interactions also becomes interesting 
from biological viewpoints. 
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